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Abstract 

One of the problems concerning entanglement witnesses (EWs) is the construction of 
them by a given set of operators. Here several multi-qubit EWs called stabilizer EWs 
are constructed by using the stabilizer operators of some given multi-qubit states such as 
GHZ, cluster and exceptional states. The general approach to manipulate the multi-qubit 
stabilizer EWs by exact (approximate) linear programming (LP) method is described 
and it is shown that the Clifford group play a crucial role in finding the hyper-planes 
encircling the feasible region. The optimality, decomposability and non-decomposability 
of constructed stabilizer EWs are discussed. 

Keywords: Entanglement Witness, Stabilizer group, Clifford Group, Linear 
Programming, Feasible Region. 
PACs Index: 03.65.Ud 
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1 Introduction 

The problem of characterizing n-qubit entanglement has motivated considerable interest in the 
literature. This problem was raised within the context of quantum information and quantum 
computation processes such as teleportation, dense coding and quantum key distribution [H [21 
[3] which consider the physical phenomenon of entanglement as a resource. Though there are 
a number of very useful and spectacular results for detecting the presence of entanglement in 
pure and mixed states of multipartite systems, the subject is still at its infancy [H El EJ [7] . 

Among the different criteria to analyze the separability of quantum states the entanglement 
witnesses (EWs) are of special interest since it has been proved that for any entangled state 
there exists at least one EW detecting it [HI [9]. The EWs are Hermitian operators which have 
non-negative expectation values over all separable states and detect some entangled states. A 
great deal of investigation has been devoted to the study of EWs, considering their decompos- 
ability, optimality [FT], optimal setups for local measurements of witnesses [T2~l [13"] and even 
their use in the characterization of entanglement in important physical systems [HI EEl [16] . 
Inside the several problems concerning the EWs, the problem of how to construct EWs by a 
given set of operators has a great importance. From a different point of view, a very useful 
approach to construct EWs is the linear programming (LP) [T71 HH [191 1201 [ZD, [22], a special 
case of convex optimization which can be solved by using very efficient algorithms such as the 
simplex and interior-point methods ( see e.g. [23, 24]). In fact, in order to a hermitian oper- 
ator W be an EW, it must posses at least one negative eigenvalue and the expectation value 
of W over any separable state must be non-negative. Therefore, for determination of EWs, 
one needs to determine the minimum value of this expectation value over the feasible region 
(the minimum value must be non-negative) and hence the problem reduces to an optimization 
over the convex set of feasible region. For example, in [T9], [20] the manipulation of generic 
Bell-states diagonal EWs has been reduced to such an optimization problem. It has been 
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shown that, if the feasible region for this optimization problem constructs a polygon by itself, 
the corresponding boundary points of the convex hull will minimize exactly the optimization 
problem. This problem is called linear programming (LP) and the simplex method is the eas- 
iest way of solving it. If the feasible region is not a polygon, with the help of tangent planes 
in this region at points which are determined either analytically or numerically, one can define 
a new convex hull which is a polygon and has encircled the feasible region. The points on the 
boundary of the polygon can approximately determine the minimum value of the optimization 
problem. Thus the approximated value is obtained via LP. In general, it is difficult to find 
this region and solve the corresponding optimization problem; thus, it is difficult to find any 
generic multipartite EW. Recently, in Ref. [21 J, a new class of EWs called reduction type EWs 
has been introduced for which the feasible regions turn out to be convex polygons. Also, in 
Ref. [22] , some kinds of Bell-states diagonal relativistic multispinor EWs have been constructed 
which can be manipulated by using exact and approximate LP method. 

On the other hand, stabilizer formalism and Clifford group operations have been proved 
to be useful in quantum error correction (theory of stabilizer codes) [251 12S1 [271 12H] , quantum 
computing, entanglement distillation [291 ESS El] and entanglement detection [T2], [13], [14]. In 
this paper, we link stabilizer theory and Clifford group operations with structure of new type 
EWs, the so-called stabilizer EWs (SEWs). As we will show all vertex points and hyper-planes 
surrounding feasible regions (i.e., the regions coming from the positivity of EWs with separable 
states) can be obtained just from a few ones by applying the Clifford group operations. The 
optimality of SEWs corresponding to hyper-planes surrounding feasible region is discussed in 
detail and it is shown that the optimality has a close connection with the common eigenvectors 
of stabilizer operators. 

The paper is organized as follows: In Section 2, we review the basic notions and definitions 
of EWs relevant to our study and describe a general approach of constructing stabilizer EWs 
by exact and approximate LP method. In Section 3, we consider the construction of SEWs 
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that can be solved by exact LP method and as instances of such SEWs, we describe the SEWs 
of GHZ and cluster states in details and give a brief discussion about SEWs of five, seven, 
eight and nine qubit stabilizer states. Also the role of Clifford group operations is studied in 
this construction. Section 4 is devoted to an analysis of optimality of the introduced SEWs. It 
is proved that some of the SEWs which correspond to surrounding half-planes of the feasible 
regions are optimal. In Section 5, we consider the decomposability or non-decomposability 
of GHZ and cluster states SEWs and show that the three-qubit SEWs are all decomposable 
but for more than three-qubit, there exist non-decomposable SEWs as well. In Section 6, we 
give some entangled mixed states that can be detected by the SEWs. Section 7 is devoted 
to construct SEWs that their feasible regions are not polygons by themselves but can be 
approximated by polygons and then solved by LP method. The paper is ended with a brief 
conclusion and two appendices. 

2 Stabilizer EWs and LP method 

2.1 Entanglement witnesses 

First let us recall the definition of entanglement and separability. An n-partite quantum mixed 
state p G B(TC) (the Hilbert space of bounded operators acting on TC = <8> ... <8> 7~td n ) is 
called fully separable if it can be written as a convex combination of pure product states, that 
is 

p = E^X^I ® k (2) X*fl ® ••■ ® l«i n) X«h (2-1) 

i 

where \cq ) with j = l,...,n are arbitrary but normalized vectors lying in the Ti^, and 
Pi > with YltiPi = 1- When this is not the case, p is called entangled. Although the 
definitions of separable and entangled states were extended to consider various partitions of 
the original system into subsystems j33j El] , throughout the paper by separability we mean 
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fully separability. 

An entanglement witness W is a Hermitian operator such that Tr(Wp s ) > for all separa- 
ble states p s and there exists at least one entangled state p e which can be detected by W, that 
is Tr(Wp e ) < 0. Note that in the aforementioned definition of EWs, we are not worry about 
the kind of entanglement of the quantum state and we are rather looking for EWs which have 
non-negative expectation values over all separable states despite the fact that they have some 
negative eigenvalues. The existence of an EW for any entangled state is a direct consequence 
of Hahn-Banach theorem [32] and the fact that the set of separable density operators is convex 
and closed. 

Based on the notion of partial transposition, the EWs are classified into two classes: de- 
composable (d-EW) and non-decomposable (nd-EW). An EW W is called decomposable if 
there exist positive operators V, Qk such that 

W = V + qT k ( 2 - 2 ) 

KcJV 

where N : = {1, 2, 3, n} and T K denotes the partial transposition with respect to partite 
K C M and it is non-decomposable if it can not be written in this form [17]. Clearly d-EW 
can not detect bound entangled states (entangled states with positive partial transpose (PPT) 
with respect to all subsystems) whereas there are some bound entangled states which can be 
detected by an nd-EW. 

Usually one is interested in finding EWs W which detect entangled states in an optimal 
way in the sense that when we subtract any positive operator from W, then it does not remain 
an EW anymore [5]. In other words, if there exist e > and a positive operator V such that 
W' = W — eP is again an EW, then we conclude that W is not optimal and otherwise it is an 
optimal EW. 
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2.2 Manipulation of EWs by exact and approximate LP method 

This subsection is devoted to describe linear programming (LP) and general approach to 
manipulate the so-called stabilizer EWs by exact or approximate LP method [23] . 
Consider a non-positive Hermitian operator of the form 



where Q t are Hermitian operators and a.'s are real parameters with a > 0. In this work, 
the operators Q i will be considered as operations of a given multi-qubit stabilizer group. The 
stabilizer operations are mutually commuting and their eigenvalues are +1 and -1. We will 
attempt to choose the real parameters Oj such that W becomes an EW. To this aim, we 
introduce the maps 



for any separable state p s . The maps P. map the convex set of separable states into a region 
which will be named feasible region. Since — 1 < P t < 1 for all i, the feasible region is bounded 
and lies inside the hypercube defined by — 1 < P i < 1 for all i. The first property of an EW is 
that its expectation value over any separable state is non- negative, i.e., the condition 



is satisfied for any point of the feasible region. For satisfying this condition, it is sufficient 
that the minimum value of J 7 W be non-negative. Therefore, for determination of EWs of type 
f)2.3p . one needs to determine the minimum value of o-o + Y^i=i over ^ ne feasible region (the 
minimum value must be non-negative) and hence the problem reduces to the optimization of 
the linear function a$ + Ym=i a iPi over the convex set of feasible region. 

We note that, the quantity T w achieves its minimum value for pure product states, since 
every separable mixed state p s can be written as a convex combination of pure product states, 




(2.3) 





^ vv :=Tr(W Ps ) = a o + ^a i P>0 
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say p s = ^jPi|Tj)(Ti| with Pi > and YliPi = 1> whence we have 

Tr(W Ps ) = ^^(WlTiXTil) > C min , (2.5) 
with C min = min 7V(W|T)(T|) 

l T >e*Vod. 

where, T) wo d. denotes the set of pure product states. In this work, we are interested in the EWs 
that their feasible regions are of simplex (or at most convex polygon) types. The manipulation 
of these EWs amounts to 

minimize F w = a + £\ a.P i ^ 
subject to E^.R-d.) >0 j = 1,2,... 

where c i . and d. are parameters of hyper-planes surrounding the feasible regions. So the 
problem reduces to a LP problem. On the basis of LP method, minimum of an objective 
function JF W always occurs at the vertices of bounded feasible region. Therefore the vertices 
of feasible region come from pure product states. 

It is necessary to distinguish between two cases: (a) exactly soluble, and (b) approximately 
soluble EWs. In the case a, the boundaries (constraints on P.) come from finite vertices arising 
from pure product states and construct a convex polygon, while in the case b the feasible region 
is not a polygon and the boundaries may be bounded convex hypersurfaces. In this case, with 
the help of tangent planes in this region at points which are determined either analytically 
or numerically, one can define a new convex hull which is a polygon encircling the feasible 
region, i.e., we approximate the boundaries with hyper-planes and clearly some vertices do not 
arise from pure product states. The points on the boundary of the polygon can approximately 
determine the minimum value of J r w in (12. 6p . Thus the approximated value is obtained via LP. 
The both cases can be solved by the well-known simplex method. The simplex algorithm is a 
common algorithm used to solve an optimization problem with a polytope feasible region, such 
as a linear programming problem. It is an improvement over the algorithm to test all feasible 
solution of the convex feasible region and then choose the optimal feasible solution. It does this 
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by moving from one vertex to an adjacent vertex, such that the objective function is improved. 
This algorithm still guarantees that the optimal point will be discovered. In addition, only in 
the worst case scenario will all vertices be tested. Here, considering the scope of this paper, 
a complete treatment of the simplex algorithm is unnecessary; for a more complete treatment 
please refer to any LP text such as [23], [21] . 

3 Exactly soluble stabilizer EWs 

In this section we consider the construction of stabilizer EWs (SEWs) which can be solved 
exactly by the LP method. In motivating this construction, we begin with EWs which can be 
constructed by the stabilizer operations of the multi-qubit GHZ state. 

But before proceeding, it should be noticed that the Hermitian operator of the form ( 12. 3j) 
can not be a SEW when all the Q/s form pairwise locally commuting set. Two operators 



are called locally commuting if [Li,Ki] = , for all i — 1, 2, n. To prove this assertion, 
consider the following operator 



Q = Li <g> ... <g> L n and Q' = K x <g> ... <8 K, 



W = aI + bQ + cQ'. 



Because of the commutativity of K± and Lj we have the 



which in turn imply that the operator W can be written as 



n 



n 



W = al + b <g)5><?|*£ ) )«'£ ) l+c 



1=1 Vi 



i=l Vi 



n 
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Now if we want W to be an EW then it must has non-negative expectation values with all 
pure product states which means that all eigenvalues (a + bX^ ...X^ + c^v} ■ ■■fJ>i$) are non- 
negative, hence W is a positive operator. Therefore, the SEWs can be constructed from the 
set of stabilizer operators Q i that at least one pair of them is not locally commuting. 

Throughout the paper, the generators of stabilizer groups are chosen according to the 
table of appendix I. Of course, this choice is arbitrary and one can take other elements as 
generators. By the method presented here we can construct SEWs (exactly or approximately) 
for completely different stabilizer groups. 

3.1 GHZ stabilizer EWs 

We consider even case of GHZ SEWs which lies in realm of exactly soluble LP problems. The 
odd case is discussed in appendix III. A similar construction can be made based on other 
elements of the GHZ stabilizer group. 

3.1.1 Even case 

Let us consider a situation in which the Hermitian operator is composed of all generators of 
GHZ stabilizer group together with all even terms 5'( ghz ) > S'_(ghz) name even re f er ^ Q ^he 
index 2k) as follows 

n n' 
k=l k=l 

where, Sj^ HZ ^ for k — 1, n are given in the table of the Appendix I and the reader is referred 
to that appendix for an overview of the stabilizer formalism. Due to the commutativity of all 
GHZ stabilizer generators, it is easy to see that the eigenvalues of W^ z are 

n n' 

«o + E(- 1 ) lfca fc + E(- 1 )' 1+i2fca ^ > V(ii > < 2> ...,i B )G{0 > l} B . (3.8) 

k=l k=l 

Evidently, when all eigenvalues are positive the above operator is positive; otherwise it may 
be a SEW. 



(3.7) 
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For a separable state p s , the positivity of 

Tr(W^ zPs ) > 
implies the positivity of the objective function 

n n' 

^GHZ = ^ + S °* P * + Yl a i,2k P l,2k ^ > ( 3 - 9 ) 
fe=l fe=l 

where 

F fc = Tr(S^ Ps ) , P h2k =Tr(SrVS^ Ps ), 

and all of the F fc 's and P 1 2fc 's lie in the interval [—1,1]. Furthermore, the operator W^ z must 
has at least one negative eigenvalue to become a SEW. To reduce the problem to a LP one and 
to determine the feasible region, we require to know the vertices, namely the extreme points 
of the feasible region. Vertex points of the feasible region come from pure product states. The 
coordinates of vertex points can take one of three values +1, -1 and 0. Regarding the above 
considerations, the product vectors and the vertex points of the feasible region coming from 
them are listed in table 1, where 



Product state 


(P2,P3,.. 


■j P'n—li P 


u Pi, Pl,2, Pl,A 


Pl,2n'-2, Pl,2n') 






(0,0,... 


,0,0, ±1,0,0, 


...,0,0) 






(0,0,... 


,0,0,0, ±1,0, 


...,0,0) 


AJ^> 




(0,0,... 


,0,0,0,0,0,... 


,0,±1) 


— 12, — / 


(_l)i a 


+ \ ...,(- 


l)^-2+^-l 5 (_ 


1)^+^,0, 0,0,. ..,0,0) 



Table 1: The product vectors and coordinates of vertices for W^ z ■ 

|*±> = |x±) 1 |x+> a |x+> 8 ...|x+> B 
A fc = ( M (2fc-i))t M(2fc) fc = 1>2 ,...,n' (3-10) 

3*,...*, = (c7< 2 >)V..(aW)*» <S>U HU) . V (<2.<3, ...,»») G {O,!}- 1 
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and \x ± ) are eigenvectors of o x with eigenvalues ±1. Here and are the phase-shift 
operator and Hadamard transform acting on particle k respectively (see appendix I). One can 
easily check by direct calculation that the convex hull of the points listed in table 1 is contained 
in the feasible region and form a (n — l)2 n ' +2 -simplex with the following boundary hyper-planes 

n' 

|P : ± P. +^(-1)^,1 = 1, J = 2,...,n, V(i 1 ,i 2l ..,i„,)G{0,ir'. 

k=l 

On the other hand, in appendix II it is shown that the feasible region is also contained in this 

simplex, i.e., the feasible region is exactly determined by the intersection of the half-spaces 

ri 

Ip^p + ^-iy^j^i. (3.ii) 

k=l 

In fact the half-spaces ( 13.111) come from the positivity of the expectation values of the operators 
2 1 3 ^ k ~ l{ h2k , J = 2,...,n , V(M 2 ,...,v)e{0,l} n 

i 2n - ± ^f Hz ) - Eti(-i)^sr z) 

over pure product states. We note that it is not necessary to consider all the above operators 
since one can obtain them just by applying some elements of the Clifford group (see Appendix 
I) on the 2n + n" (compare with (n — 1)2™' +2 ) following operators 

I 2n ± (SfOHZ, + 5 (GHZ) + sJGHZ)) j = , ^ ^ 

I 2n - S[ G ^ - Sgf> - Eti 4T Z) 3 = h n. 
For example we get the operator S = I 2n + S[ cnz) - S^ GRZ) - S[ G f z) + Ylk+j S f^ Z) from 
the operator S' = I 2 „ + S[ GHZ) + S^ GHZ) + YJk=i S f^ Z) under conjugation with the Clifford 
operation i.e., 

Now the problem of finding a pre-SEW (a hermitian operator with non-negative expectation 
value over any separable state) of the form (13. 7p is reduced to the LP problem 

minimize = a + £Li a k P k + J2t=i a h2k P i, 2k 

subject to \P 1 ±P j +Y;ti(-l) ik Pi,J <1, J = 2,...,n, V (n, z 2 , v) 6 {0, 1}-' 

(3.13) 
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On the basis of LP method, minimum of an objective function always occurs at the vertices 
of the bounded feasible region. Hence, if we put the coordinates of the vertices (see table 1) 
in the objective function (13. 9p and require the non-negativity of the objective function on all 
vertices, we get the conditions 

a > , a >KI , «o>Er= 2 kl 

(3.14) 

%>\a 12k \ k = l,...,n' 
for parameters dj. Evidently, these conditions are sufficient to ensure that the objective func- 
tion is non- negative on the whole of the feasible region. If we take a = (n — 1), a k = — 1, 
for all k = l,...,n and a 12k = 0, for all k = l,...,n', which fulfill all the conditions of Eq. 
(1513)1 . then we get the SEW stated in Eq. (21) of Ref. PE|. Also by taking a = 1, a, = -1 
and a m = a lm = —1 (m > 2 is even) we have the SEWs stated in Eq. (21) of the mentioned 
reference as special cases. 

Fixing a in the space of parameters, all of the czj's lie inside the polygon defined by 
inequalities (13.141) . Now in order that the operator of Eq. (13.71) becomes non-negative, all of 
its eigenvalues in ( 13.8f) must be non-negative. The intersection of half-spaces arising from the 
non-negativity of the eigenvalues form a polyhedron inside the aforementioned polygon. The 
complement of this polyhedron in the polygon is the where that the operator (13. 7p is SEW 
and will be named the SEWs region. 

We assert that the SEWs region is non-empty. To confirm this assertion, we discuss the 
case that all parameters a, are positive since the discussion for other cases can be easily come 
from by replacing any parameter by its negative value (except a which is always positive). 
Because of the symmetry between the parameters a 1 and a 2k 's (A; = 1, ...,n'), we can assume 
without loss of generality that a 2 > a 4 > a 6 > ... > a 2n ,. With this assumption, all of the 
2 n " + n' eigenvalues (with n" = f 1 ^]): 

% + a, + Y%Li a 2j + YTjUi-^^^+x + Elti a h 2 k V («3, is: W'+i) e {0, l} n " 
a + a 1 — a 2l + ^^ fc=1 a 12k + z\2j=2 a 3 I = I, ...,n 
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are non-negative and each of the 2 n — (2 n " + n') remaining ones can take negative values. 
For example, consider the Hermitian operator 

W<2, = a J 4 + a.Sf^ + a 2 Sf HZ > + a U2 S^Sf HZ > (3.15) 

with the following eigenvalues 



(3.16) 



u 1 = a + a, + a 2 + a lj2 , cu 2 = a + a, - a 2 - a 12 
a; 3 = a - a, + a 2 - a 12 , w 4 = a - a, - a 2 + a 12 . 
We need only to consider the product state \x + ) \x + ) corresponding to the vertex point (1, 0, 0) 
since the product states corresponding to the other vertex points can be obtained by applying 
the Clifford operations H ® H , M <g> M and o z ® I on this product state. Putting the vertex 
points in Tr(W^ z |T)(T|) > yields 

% > Kl, a > \a 2 \, a > \a 12 \ . 

So in the parameters space, the allowed values of a's lie inside a cube with edge length a . The 
intersection of half-spaces uj i > (i — 1, ..,4) is a polyhedron inside the cube whose vertices 
coincide with four vertices of the cube and contains just the positive operators; the remaining 
part of the cube is the region of SEWs. On the other hand the variables P i lie in the interval 
[—1,1] and form a cube in the space of variables. The convex hull of vertex points lies inside 
this cube and has the eight boundary half-spaces 

\P,±P 2 + P lt2 \ < 1 , \P X ±P 2 -P lt2 \<\. (3.17) 

The above half-spaces define the feasible region (see Fig.l). Four of these half-spaces which 
correspond to the positive operators 

l V GHZ = (h + + 5 (GHZ) + sgmz)) = 4 

2 V GHZ = ^CV GHZ )a[ 1] =h~ ^ (GHZ) + *f HZ) - ^ (GHZ) 

3 -p ghz = °^v GH M? = h + ^ (GHZ) - * 2 (GHZ) - ^ GHZ) 

A V GHZ = <7UCV GHZ )aM = J 4 _ 5 (G HZ) _ S (GHZ) + 5 (GHZ) 
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are in one-one correspondence with four vertices of the cube in parameter space which are the 
same as the vertices of polyhedron formed by the positive operators. 
For the purpose of later use, we introduce 

= w* 1 ® (^r ® ... ® wivw..o>. (3.i9) 

where |Vv..o) = ^=(|00...0) + is the n-qubit GHZ state. As implied by the Eq. (Engl) , 

the three last positive operators can be obtained from the first one via the action of some 
operations of the Clifford group. The other four boundary half-spaces which correspond to the 
optimal d-EWs 



X W£S = h- Sf HZ ) - HZ ) - S( 2 GHZ ) = 4(1^)^ If 1 

W£*5 = CrW^W^ )0"W = J 4 - 5(GHZ) + ^(GHZ) + £-(GHZ) 



(3.20) 

3yy(opt) _ a (l)nyy(opt)\ a {l) = J + ^(GHZ) _ ^(GHZ) + ^(GHZ) 

GHZ z * GHZ ' z 4 1 2 12 

4yyfcp*) = o-W^W^^crW = 7 4 + S'f 0112 ) + S'p 112 ) - S'(ghz) 

are in one-one correspondence with the remaining four vertices of the cube in parameters space. 
From Eq. (13.201) we see that the three last optimal d-EWs can be also obtained from the first 
one via the action of some operations of the Clifford group. So as we had in [TU], the operators 
corresponding to the boundary planes are either optimal SEWs or positive operators. In this 
case, all of the witnesses are d-EWs since we can write them as a convex combination of an 
optimal d-EW and a positive operator from its opposite positive boundary plane. 

3.2 Multi-qubit cluster EWs 

We continue with EWs which can be constructed by the stabilizer operators of the cluster 
state and again consider two even and odd cases of the cluster SEWs which lie in the realm of 
exact LP problems (refer to appendix III for odd case). 
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3.2.1 Even case 



n + 1 



1 (3.21) 



Let us consider the following Hermitian operators 

n' 

= + E ^ + a ^S£l + a 2m _ lt2m S£lS^, m = 2, 

k=l 

In addition to the above operators, one can consider other Hermitian operators which dif- 
fer from the above operators only in the last terms, that is the last terms of them are 
a, ,Sp\Sp\ with m = 2,..., T^^-l. However, we will consider only the operators 

2m — 2,2m- 1 2m-2 2m— 1 ' ' L 2 J ' J ^ 



( I3.2ip since the treatment is the same for others. Due to the commutativity of all cluster 
stabilizer generators, it is easy to see that the eigenvalues of W^ 1 ' are 



2m — 1,2m 



, V(M 2 ,...,z n ) G {0,1}" (3.22) 



a o+E(- 1 ) i2Ja ^+(- 1 ) i2m " la -- 1 +(- 1 ) i2m " 1+i2m 

3=1 

To reduce the problem to a LP one and determine the feasible region, we require to know 



the vertices, namely the extreme points of the feasible region. For a separable state p s , the 
non-negativity of 

Tr(W^p s ) > 
implies the non-negativity of the objective function 

~n + 1" 



^2m-l,2m^2m-l,2m ' m 2,..., 



1 (3.23) 



J~ — a + a 2fc^2fc + a 2m-1^2m-l + 

k=l 

where, 

P 2k = Tr(S^ Ps ) , F 2m _ li2m = Tr(S^_A^Ps), 
and all of the -P 2fe 's and -P 2m _ 12m 's lie in the interval [—1,1]. The product vectors and the vertex 
points of the feasible region coming from these product vectors are listed in table 2 
where 

|$) = \z + ) 1 \x+) 2 \z+) 3 \x+) 4 \z+) 5 ...\x+) n _ 1 \z + ) n 

, , = ®U , V (i u i 2 , v) G {0, 1}"' 

1 ' ' ' ' ' TV J 



A'M = ff(2m-2) ff(2m-l) jj(2m) 

ii / £l .... ,2 / 

1 ' ' 77/ 1 ' ' TV 



A 



y^(ei)) jj{2m-2) ^{2m-l) jj(2m-l) j^(2m) 
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Product state 


(-P2) -^4, • * • 5 -P2m-4, -P2m-2, -P2m-lj ^mi -P2m+2i •••) -^Wi', -P2m-l,2m) 


AM 1$) 


((_!)* (_i)^- 2 , (-1)^,0, (-1)*™, (-1)^+S (-1)V,0) 


il,-»V 


d>) 


((_!)« (_i)« (-1)^,0, ±1,0, (-1)W S (_i)v 0) 


n.-.V 




((-1Y\ ( — l)* 2 , (-1)*- 2 , 0, 0, 0, (-1)*™+!, (_1)V, (_l)*m) 



Table 2: The product vectors and coordinates of vertices for 

For a given m, the convex hull of the above vertices, the feasible region, is a (2n' + 12)-simplex 
formed by the intersection of the following half-spaces 

\P2m-l i P2m-2 + P2m-l,2m\ < 1 
lAm-l ± P2m-2 ~ P2m-l,2m\ < 1 

|-p2m-l ± -P2m + P2m-l,2m\ < 1 (3.24) 

l-fW-l ± -P2m — -P2m-l,2m| < 1 

\p2k\ ^ 1 j m,m — 1 ^ k — 1, ...,n' 

(see Appendix II). In fact the half-spaces (I3.24p come from the non- negativity of the expectation 
values of their corresponding operators 

I + Sf > ± SP , + SP , Sf) , / - Sf) =f Sf) , - Sf) 

2m— 1 2m — 2 2m — 1 2m ' 2m— 1 1 2m — 2 2m— 1 2m 

/ + Sf > ± SP - SP SP , / - SP , t SP , + SP SP 

2m— 1 2m — 2 2m— 1 2m ' 2m-l 1 2m — 2 2m— 1 2m 

i + sp ±sp + sp sp , i-sp T sp-sp sp 

2m — 1 2m 2m— 1 2m ' 2m — 1 2 m 2m — 1 2 m 

i + # c > ±sp-spsf> , /-sp Tsp + spsp 

2m— 1 2m 2m— 1 2m ' 2m— 1 1 2m 2m— 1 2m 

i" ± S^ ) , m, m — 1 ^ k = 1, .... n' 

2fc ' ' ' ' ' 

over pure product states. We note that it is not necessary to consider all the above operators, 
since one can obtain them just by applying some elements of the Clifford group on the 4 
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(compare with 2n' + 12) following operators 



i±sp ±sp±sp sp 

2m — 1 2m 2m— 1 2m 

/_5(C) 5jc) S (C) s p (3.25) 

2m — 1 2m — 2 2m — 1 2m 

For instance, the Clifford operation 

U = (CN i2 )(CN 53 )(CN 13 )(CN 2i ) e Cl(n) 
transforms Sp to Sp by conjugation, i.e., 

usfW = sp 

2 4 

Now the problem of finding a pre-SEW of the form (I3.2ip is reduced to a LP problem with 
objective function (13.231) and constraints (13.241) . If we put the coordinates of vertices (see table 
2) in the objective function (I3.23f) and require the non-negativity of the objective function on 
all vertices we get the conditions 

% > Ejli K-l 

a <> > E7=i 2 \ a 2j I + Ejlm+i \ a 2j I + k m _! I ( 3 - 26 ) 

^ > \--\m—2 I I , v-^ n ' | lil I 

a — 2—lj=\ \ a 2j \ + 2-lj=m+l \ a 2j \ + l a 2m-l,2ml 

for the parameters Oj. Evidently, these conditions are sufficient to ensure that the objective 
function is non-negative on the whole of the feasible region. Cluster SEWs (I3.2ip and the odd 
case discussed in appendix III contain the SEWs in Eqs. (36) and (37) of Ref. [10] as special 
cases. 

Fixing a in the space of parameters, all of the a's lie inside the polygon defined by inequal- 
ities ( 13.26H . Now in order that the operator of Eq. ( 13.211) becomes positive, all of its eigenvalues 
in (13.22p must be non-negative. The intersection of half-spaces arising from the non-negativity 
of eigenvalues form a polyhedron inside the aforementioned polygon. The same reasoning as in 
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the even case of GHZ SEWs, shows that the SEWs region is non-empty. For example, consider 
the operator 

= a / 4 + 0i SP + a 2 SP + a 4 Sp + a^Sfl 
The eigenvalues of this operator are 

lo 1 = a + a 1 +a 2 + a A + a 12 , uj 2 

u 3 — a + a x + a 2 — a 4 + a lj2 , uj a 

u 5 = a - a x + a 2 + a 4 - a 12 , u 6 

u> r = a — a x — a 2 + a 4 + a 12 , a; g 

Without loss of generality we can assume that a x > a 2 . With this assumption, the first 
four eigenvalues u 1 ,u 2 ,uj 3 and u 4 are always positive. Now let u 5 and a; 6 be negative, i.e., 
a + a 2 < a x + a 12 . In this case, u 7 and a; g can not be negative and vice versa. Therefore 
with these considerations, among the eight eigenvalues only the pair u 5 ,u 6 or coV,w 8 can be 
negative. The explicit form of some four-qubit cluster SEWs is postponed to section [51 

4 Optimality of SEWs 

Another advantage of stabilizer EWs is that the optimality of the EWs corresponding to 
the boundary hypereplanes of feasible region can be easily determined by a simple method 
presented here. Consider an EW corresponding to one of the hyper-planes in which three 
terms S v S j and S.S. appear simultaneously such as 

w = / + (-i)* 1 ^ + (-irs;. + (-ly^s; + ... v i u i 2 , i 3 , ... e {o, i}. (4.27) 

If there exist e > and a positive operator V = |?/>) (■!/>!, such that W' = W — e\i/j) is again an 
EW then we conclude that W is not optimal, otherwise it is. Note that there is no restriction 
in taking V as a pure positive operator since every positive operator can be expressed as a sum 
of pure positive operators with positive coefficients, i.e., V = AJ^i)^! with all > 0. If 



= a + a 1 — a 2 + a 4 — a 
= a + a 1 — a 2 — a 4 — a 
— a — a 1 + a 2 — a 4 — a 
= a — a x — a, 2 — a 4 + a 
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W" be an EW, then \ip) has to satisfy the constraint Tr( |^)(^|Y)(T|) = | | TC) | 2 = for any 
pure product state |T) satisfying Tr(W|T)(T|) = 0. In other words, has to be orthogonal 
to all such pure product states. 

Since in SEWs of the form (I4.27P considered so far there is no pair of locally commuting 
operators, it is always possible to find pure product vectors |T) for which one of the relations 

SJT) = (-l^T) 

S.\T) = (-l) i2+1 |T) (4.28) 
^|T) = (-1)^ +1 |T) 

hold. The expectation value of W over such |T)'s is zero and \ip) cannot contain such pure 
product vectors. All the eigenvectors of a stabilizer operation can be chosen as pure product 
vectors, half of them with eigenvalue +1 and the other half with eigenvalue -1, such that the 
expectation value of other stabilizer operations over them be zero. Because of Hermiticity of 
stabilizer operations, their eigenvectors can be used as a basis. 

Let us assume that |Tfc)'s are pure product eigenvectors of S i with eigenvalues (— l) n+1 
and |T^)'s are its pure product eigenvectors with eigenvalues (— l) n that have been chosen 
according to the above prescription. So the expectation value of W over |Tj.)'s is zero and 
cannot contain |Tfe)'s that is = J2k\"^k)- This implies that S^ip) = (— By 
the same reasoning we conclude that SAi/j) = (— lY 2 \ip) and S^Aip) = (— On the 
other hand, we have S^^i/j) = (— l)* 2 ^!^) = (— l) n+ * 2 |^). Hence, if ^ %\ + 12, i.e., if 
i 3 = i\ + %i + 1, we get into a contradiction and W is optimal. Therefore, among all SEWs of 
the form (14.271) the following ones are optimal 

yv op t = 1 + (-i)* 1 ^ + (-iy 2 s j + (-iy^ +1 s i s j + ... v (i 1} i 2 , ...) e {o, (4.29) 

With the same reasoning as above one can conclude that any SEW of the general form 

W = 1 + (-IT'S, + (-l) 42 ^. + (-iy 3 SA + ... V h, i 2 , ... e {0, 1}. (4.30) 

with j ^ k is not optimal. 
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For instance, in the case of the three-qubit Gif Z-state, 

WSL = a oI 8 + «^; GHZ) + a 2 Si G HZ ) + a 3 Sf HZ ) + a^Sf^Sf^ (4.31) 
the boundary half-spaces of the feasible region are 

l^iP.+Pj <1 , l^iP.-Pj^l j =2,3 (4.32) 

Using Clifford group operations, we can obtain all of these half-spaces only from the three 
half-spaces 

\P : +P 2 + PJ <1 , P 1 + P 3 + P 1 , 2 <1 (4-33) 
The operators corresponding to the above boundary half-spaces are 

Q GHZ = h + sr z) + * 2 (GHZ) + s( GHZ >s( GHZ > = 4(|^ 000 )(^ 00 | + IVvX^U) 

2 W CTZ = J 8 - 5f HZ ) - 5| GHZ ) - ^hz )5 (ghz) = 4|^ loi )^ i01 | +4(|^ 110 )(^ 110 |) T2 . 
It is seen that, in agreement with the above argument, 1 W G/J - Z is an optimal SEW but 2 W GHZ 
is not. Also, for the three-qubit cluster state, 

W'f = o / 8 + a.SP + a 2 SP + a 3 Sp + a^S^S™ (4.35) 

using Clifford group operations, we can obtain all of the boundary half-spaces only from the 
three half-spaces 

\Pi + P 2 + Pi*\ < 1 , P 2 + P 3 + P^<^ 
and the operators corresponding to the above boundary half-spaces are 

hWhwq ghz hWhw = i s + sp + sp + spsp 

HWhWW ghz HV>H<® = I 8 - Sp - Sp - SPSP (4.36) 

hWhv> 2 w ghz hWhv> = i s - sp - sp - spsp 

Clearly, local unitary operations Ui oca i do not change the optimality of EWs under the conju- 
gation action such as UiocaiW op U} ocal , hence among the above operators, the second one remains 
optimal while the third one remains non-optimal. 
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5 Decomposability of SEWs 

Another interesting feature of EWs which is necessary to study about SEWs is decomposability. 
As it is well-known that every two-qubit EW is decomposable [HI El US], we discuss the three- 
qubit systems or more. 

5.1 Decomposability of 

First consider three-qubit GHZ SEWs. The inequalities (13.141) show that in the space of 
parameters all GHZ SEWs lie inside the hypercube (again by fixing a ) but this statement 
does not mean that any point of the region inside the hypercube is an SEW. The region defined 
by the inequalities 

a + (-l) 11 ^ + (-l) i2 a 2 + (-l) i3 a 3 + (-l) il+ ^a lj2 > (h, i 2 , z 3 ) G {0, l} 3 (5.37) 

is the place inside the hypercube where the operator have just positive eigenvalues and 

hence is positive. First we consider the decomposability or non-decomposability of SEWs 
lying over the edges of the hypercube. These SEWs come from 1 W GHZ and 2 W GHZ of ( 14. 34ft 
by Clifford operations. The l V^ GHZ and SEWs coming from it are optimal decomposable since 
their partial transpositions with respect to some particles are positive. 

Now in the space of parameters a, we consider the coordinates of points 
Putting the following four points (which lie over the edges of hypercube) in W G 3 2 Z gives the 
following optimal SEWs 

(1, 1, 0, -1) - I 8 + S[ GH ^ + HZ > - SJGHZ)£(GHZ) 
(-1,1,0,1) - I g — ^(GHZ) _|_ g(GHZ) _|_ g(GHZ)g(GHZ) 
(1,-1,0,1) -> / 8 + ^ GHZ ) - Sf HZ ) + iS (GHZ) 5 (GHZ) 
(_1 ; -I, 0, 1) -> I 8 - S[ GU ^ _ S-(GHZ) _ 5 (GHZ) S (GHZ) 

The above SEWs are optimal decomposable since their partial transpositions with respect to 
some particles are positive . A convex cone which may be formed by connecting every four 
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points of Eq. (15.381) to its opposite positive hyper-plane in Eq. (15.371) is d-SEWs. Note that 
the remaining operators in Eq. (14.321) coming from some points in the space of parameters are 
either d-SEW or positive. Therefore we conclude that all the three-qubit GHZ stabilizer EWs 
are decomposable. The discussion for more than three-qubit is rather complicated. It is clear 
that every EW with positive partial transpose with respect to some particles is decomposable. 
Therefore imposing the condition 

n 

ao+E(- 1 ) Jfc ^ + E(- 1 ) n+ifc+ X fc + E (-ir + x fc > ° (h,..,in)e{o,ir (5.39) 

k=l keB keA\B 

which in turn implies W Te > 0, yields the GHZ decomposable SEWs where the B is any 
nonempty subset of the set A = {2,4, ...,2n'}. Here taking partial transpose with respect to 
the particles 2j and 2j — 1 with j = 1, n' leads to the same result. 

In order to show that the for n > 4 contain some nd-EWs, we discuss the four-qubit 

case in detail. From ( 13. 71) . we have 

k=l k=l 

Using the local Clifford operations, all the 48 Hermitian operators corresponding to boundary 
half-spaces of the feasible region can be obtained only from the following 5 ones 

IWgk = J l6 + SJGHZ) + SfHZ) + 5 (GHZ) S (GHZ) + S (GHZ) £-(GHZ) 

2 W( 4 2 z = I 16 + S[ GHZ) + S[ GHZ) + S[ GHZ) S^ 112 ) + 5JGHZ) ^(GHZ) 

3yyg)^ _ _ £(GHZ) _ ^(GHZ) _ £(GHZ) £(GHZ) _ g(GW)g(GW) (5.40) 
4 )/y( 4 ) = I — S[ GnZ ^ — 5*( GHZ ) — g( GHZ ) g( GHZ ) _ g(GHZ) g(GHZ) 



5 



yy(4) _ j _ g(GHZ) _ g(GHZ) _ g(GHZ) g(GHZ) _ g(GHZ) g(GHZ) 

Now consider the following density matrices 

p± = 1 [J ib ± ^(Sj GHZ ) + gfHZ^GHZ) + gfHZ^GHZ) - S-JGHZJ^GHZJ^GHZ))] (5 4l) 

One can easily check that p ± are PPT entangled states and can be detected by the above 
SEWs, i.e., 

TrfW^p.) = ~ for 2 = 1,2, (5.42) 
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and 



Tr ^alzP + ) = ~\ for ^ = 3,4,5 (5.43) 



which means that all SEWs stated in Eq. (15.401) are nd-SEWs. On the other hand, by the 
(Ea , and are optimal SEWs. 

GHZ GHZ 

Moreover, by the following transformations 

— W = U local W^Ulcal 



(5.44) 

± r H± ~ w local H± w local 



P± * P' + = Ulocalpjj} 



where Ui oca i may be any local unitary Clifford operation we can get the new nd-SEWs 'W 
which can detect the PPT entangled states ft ' . It is necessary to mention that local unitary 
operations transform a PPT entangled state to a PPT one. 

5.2 Decomposability of W^ 1 ' 

Since the three-qubit cluster SEWs are transformed to three-qubit GHZ SEWs by local unitary 
Clifford operations as in Eq. (14.361) therefore they are also d-SEWs. For more than three-qubit 
the discussion is similar to the GHZ one. The SEWs 

I _ Q(C) _ C(C) _ o(C) S (C) (5 _ 45) 

2 n 2m-l 2m 2m-l 2m \ / 

are optimal d-SEWs since they have positive partial transpose with respect to the particle 
2m — 1 or 2m. Again a convex cone which may be formed by connecting every points of Eq. 
(I5.45P in the space of parameters to its opposite positive hyper-planes 

n' 

^ + E^ 1 )^ + (- 1 ) l2m_la 2m- 1 + (-ly^+^m-Mm > (5.46) 

3=1 

for all (ii, i 2 , in) G {0, l} n , are d-SEWs. 

For illustration, we discuss the odd case of 4-qubit cluster SEW in detail. From flIII-61) . we 
have 

<> = «./„ + E<W^2, + "^f ' + a„S<°>Sf > 

k=0 
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Using the local Clifford operations, all the 14 Hermitian operators corresponding to boundary 
half-spaces of the feasible region can be obtained only from the following 3 ones 

^wf = / 4 + sp + sp + spsp, 

2>v'( 4 ) = I - £(c) _ 5(c) _ 5 , (c) i5 (C) ! (5.47) 

2^ 2 3 2 3 

*wf = j — sf ) - sp - spsp. 

O 2 J- ^ ^3 

Among the above operators, 1 W'[5 > is positive since 1 W'^ = ( I + Sjp)(I + SP), and if we 
take partial transpose of the second one with respect to second particle we get 

CWff* =(I- SP)(I - SP) > 0, (5.48) 

so 2 Wf is an optimal d-SEW. 

Although we could not find bound entangled states which can be detected by exactly soluble 
cluster SEWs however we will be able to find such entangled states for approximately soluble 
cluster SEWs as discussed in section [7] and therefore we postpone to subsection 17.21 for more 
details. 



6 Separable and Entangled stabilizer states 

Once again consider the general form of operators which is the same as Eq. (12.31) . i.e., 

i 

P ■= E K^ n -A jl Sj\..St- k = o,/.,.. + £ CjAj (6.49) 

h j2,-Jn-fc = j^O 

where for simplicity we have renamed the S^S^ 2 ...S J n n _~ k and ^ 1 ,j 2 ,...,j n _ fc by Aj and Cj respec- 
tively. Positivity of p together with &o,o,...,o = Co = ™ make (16.491) a density matrix. On the 
other hand, we assert that the conditions 

1 1 1 

jlj2,-jm=0 j^O 

yields separable state. To see this we note that for any element Aj of S n -k the operator I + A. 
is separable because it is the projection operator on the space spanned by the pure product 
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eigenvectors of A j corresponding to the eigenvalues +1. So any convex combination of the 
operators I + A such as 

e~, ■= ^ + ^ E^ + A) = ^ + ^ E^i ( 6 - 51 ) 

is separable where ^2j-^ Pj = 1 and < /x < 1. The same statement holds if we replace some 
I + A. by / — A j in the above equation. Now if we consider all Cj to be positive in Eq. (I6.49P 
and rename Z^ Pj by Cj (with j ^ 0) we conclude that the condition (16.501) is satisfied and 
therefore p is separable. For the cases that some Cj are negative it is enough to replace some 
I+A. by I — A. in the Eq. (16.511) and proceed the same way as described above. Consequently 
we get a family of separable states expressed in terms of the elements of the stabilizer group 
provided that the condition (16.501) satisfies. In the following, some entangled states including 
PPT ones which can be detected by GHZ and cluster SEWs are introduced. 

6.1 Entangled states which can be detected by 

& J GHZ 

Now we assert that GHZ stabilizer EWs can detect some mixed density matrices. To this aim 
consider the following operator 

& := £ b^,...j n S[^S^...S^ (6.52) 

jl,32,-,3n=0 

which due to tracelessness of ,S( GHZ ) the condition Tr(pW^) = 1 gives 6 0j0j-- . i0 = ^ an d the 

positivity of density matrix impose 
i 

E (-l) iljl+i2h+ - +injn b hda _j n > , V(*i,i 2 ,...,i B )e{0,l} B (6.53) 

iij2,-,>=o 

to its eigenvalues. An interesting case is when all coefficients are equal to &j 1)3 - 2 ,...j n = 4r which 
is coincides with the n-qubit GHZ state 

n 

\A ..,)(Ao..,\= ¥ u( i+s r z) y 

3=1 
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This density matrix has 2 n terms which except b the other are arbitrary parameters with 
the constraints in Eq. (16.551) . These 2 n constraints forms a simplex polygon in a 2 n — 1 
dimensional space with coordinate variables bj lt j 2j _j n (excepted & 00 _ ). Furthermore if we 
want p^\ z becomes a PPT entangled state in the sense that its partial transpose is positive 
definite with respect to any particle, i.e., p^ z Tl — ^ w ^ n i = 1, w then we must have 
i 

E {' ')' + 6 0> ,w«} >o -V (/,./,. i n ) G {0,1}" 

Jl,i2,—,in=0 

Introducing the new parameters b i = b 1 with a 1 in the zth position, and b 1 . = 
6 1 1 with a 1 in the jth position , and using the orthogonality ( 11-21) of S/s, then the 
condition for detectability of p( n ' by yy( n ) can be written as 

GHZ GHZ 

n n' 

Tr(W^M z ) = ^ + E a A + E a ^h™ < 

k=l k=l 

6.2 Entangled states which can be detected by W^ n) 

Now we assert that the above cluster stabilizer EWs can detect some mixed density matrices. 
To this aim consider the following operator 

E b hAr .. dn Sp h Sp h ...S^ (6.54) 

which due to traceless of the condition Tr(p) = 1 gives b = and the positivity of 

density matrix impose 
l 

E b juj2 _ jn > , V (h, k, i n ) e {0, IT (6.55) 

Jl J2v ,jn=0 

to the its eigenvalues. An interesting case is when all coefficients are equal to bj u j 2t _j n = ^ 
which is coincides with the n-qubit cluster state 

1 n 

i c >< c i=^n( /+ ^ c) )- 
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In order to can be detected by an odd case we must have 

n 

TrQtffpW) = ^ + J> 2fc+1 6 2fc+1 + a 2m b 2m + a 2mM b 2mM < 0. 

k=i 

7 Approximate stabilizer EWs 

So far, we have considered SEWs which can be exactly solved by LP method. In this section, 
we consider approximately soluble SEWs which come from by adding some other members 
of stabilizer group to exactly soluble SEWs. In all of the SEWs discussed in section [3j the 
boundary half-spaces arise from the vertices which themselves come from pure product states 
and the resulting inequalities did not offend against the convex hull of vertices at all. But 
by adding some terms to exactly soluble SEWs, it may be happen that the feasible region be 
convex with curvature on some boundaries and the problem can not be solved by exactly LP 
method. In these cases the linear constraints no longer arise from convex hull of the vertices 
coming from pure product states. Hence we transform such problem to approximately soluble 
LP one. Our approach is to draw the hyper-planes tangent to feasible region and parallel 
to hyper-planes coming from vertices and in this way we enclose the feasible regions by such 
hyper-planes. It is clear that in this extension, the vertices no longer arise from pure product 
states. 

7.1 Approximate n-qubit GHZ SEWs 

For the even case of GHZ SEWs we add one of the statements sCGhzj^ghz) n _ //\ t 

1 21 + 1 V ' ' / 

Eq. ([32D as 

n n' 

WW = a I 2n + y a fc S< GHZ > + y a, 2 ^ GHZ ^ GHZ > + a 2l+1 S^S^ (7.56) 

GHZ (ap) ' J k / J 1 - 2k 1 2k 1,21+1 1 21 + 1 v > 

k=l k=l 
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and try to solve it by the LP method. The eigenvalues of Wi™* 1 are 

J J ^ GHZ (ap) 

n n 1 

«o + + Et- 1 )^^ + (- 1 ) tl+ ^ +la i, 2i+ i . V <2, ^n) G {0, 1}" 

j=l fc=l 

The coordinates of the vertices which arise from pure product vectors are listed in the table 3 



Product state 


\Pli P31 ■■■■> Pn-1, Pn, 


Pi, Pl,2, Pl,4, ■■ 


j Pl,2n'-2, P\,2n'i P\,2l+\) 




(0,0,. 


..,0,0, ±1,0,0, 


...,0,0,0) 




(0,0,. 


..,0,0,0, ±1,0, 


...,0,0,0) 




(0,0,. 


..,0,0,0,0,0,.. 


,0,±1,0) 


j\(ap)|\j/±\ 

21 + 1 ' 1 


(0,0,. 


..,0,0,0,0,0,.. 


,0,0,±1) 


^i 2 ,...,i„\^ + ) 


((-i)\(-i)v*s,...,(. 


_l)^-2+^-l 5 (_ 


-1)^-1+^,0, 0,0,. ..,0,0,0) 



Table 3: The product vectors which seem to be the vertices for Wl") 

GHZ (a v ) 



where 

A (ap) = (M (2Z) ) t M (2i+1) . 

Choosing any N\ = n + n' + 1 points among N 2 = 2 n_1 + n' + 2 above points we get the 
following half-spaces 

ri 

\P, ±P+ E(" 1 ) lfcp i,- + (- 1 ) V+ip 1 , 2i+1 1 < ^ =?, J = 2, n, (7.57) 
k=i 

where v+i) G {0,1}"' +1 . But calculations show that the inequalities offend against 1 

up to iimax = (see appendix II). This shows that the problem does not lie in the realm 
of exactly soluble LP problems and we have to use approximate LP. To this aim, we shift 
aforementioned hyper-planes parallel to themselves such that they reach to maximum value 
Umax = ■ Ori the other hand the maximum shifting is where the hyper-planes become 
tangent to convex region coming from pure product states and in this manner we will be able 
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to encircle the feasible region by the half-spaces 



1^ + + Eli + *Wl < ^ 3 = 1, 

I P + P + Y^ n ' P IP I < l + V^ 

I" 1 1 ' ± 21 + 1 ' Z^lk=l 1,2k ' ± 1,21+1 I — 2 

Pi + P v+1 + Eti Px.„ < ^ ^ J = 1, - n" 

P, < 1 



(7.58) 



where again we have used the Clifford group and write just the generating half-spaces. Due to 
the above inequalities the problem is reduced to the LP problem 



minimize 



rr(ww l 7 )( 7 |) 



(up) 



s.t. 



\Pi±P + Zti(-±Y k Pi,2 k + (-i^'+^+J < 



1+V2 



j = 2,...,n 



P < 1 



Pl,2 k < 1 



P i,n+i < 1 



i — 1, n 
k = 1, n' 



for all (ii, Z2, v, V+i) G {0, where it can be solved by simplex method. 
The intersections of the half-spaces in the above equation form a convex polygon whose 
vertices lie at any permutation P[,P'p P[ 2k {k = l,...,n') and P[ 2l+1 with a given j (j = 
2, ...,ri) of the points listed in table 4 where P' 's are defined by shifting the P 's for all 
(ii, ...,i n ,k,2, —,ii,2n>) e {0, 1}"+" . 

So in order that the expectation value of VVl") be non-negative over any pure product state, 



(ap) 



the following inequalities and any inequality obtained from them by permuting the parameters 
a 11 a j1 a 12k (k = 1, ...,n') and a 12l+1 with a given j for j = 2, ...,n, must be fulfilled 

% + ELi(-i) ifc « fc + Elit-i)^ 2 ^ + ^i,2 l+ i > o 

such that (-1)^ + (-!)** + Efelil- 1 )* 1 ' 2 " = 2 



« + ELi(-i) J X + EtiC-i)* 1 -"^ + 



3-Vg, 



1,21+1 



> 



such that (-1)" + + Efelil- 1 )^' 2 " = -2 
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p' p' p' p' p' p' 


p/ p/ \ 

^l,2n'' - r l,2Z+l J 






(7—1 V2 


C 1 Va 


C 1 W- 


1 j <L^ik=l i,2fc 


( 1 Vl 2n'-2 


f 1 Vl 2n' 


2 J 




( 1 


C 1 V»- 

•••5 I X J 


9P' + P' + V™' ,P' = -2 


( lVl2n'-2 

••) V 1 


f 1 Vl 2n' 


2 J 


(7— iV 2 


C— 1 V 3 


C— 1 W- 


3 P' + p' + J2 n 'p' =1 


C_1 Vl,2n'-2 


(—] Vl,2n' 


\/2-n 
2 J 


((-l) i2 , 


(-!)*», 


...,(-!)*»- 


- 1 , . 
3 p' + p' + yi I?' =-1 


„,(-l)*i,w- aj 




l-%/2\ 
2 J 



Table 4: The coordinates of vertices for 

GHZ (ap) 

« + ELi(-i)^ + EtiC-i) 41 ^ + ^«i, 2;+ i > o 

such that (-1) 41 + (-1)^ + Efe=i(- 1 ) il ' 2fc = 1 

« + ELi(-i)'*^ + Eti(-i) !l ' 2 ^i,, + ^k, 2;+ i > o 

(7.59) 

such that + (-l)v + Efc=i( -1 ) il,a * = _1 

Similarly, one could repeat this approximation for the odd case of GHZ SEWs like above. 

As the case of exactly soluble GHZ SEWs, we assert that there exist some nd-SEWs among 

the approximately soluble GHZ SEWs. To see this, consider the four-qubit GHZ SEWs 

W ± = 1 + ll J ig + 5<(GHZ) + ^(GHZ) + S (GHZ) ^(GHZ) + ^(GHZ) g(GHZ) ± g(GHZ) ^(GHZ) ^ 

which both can detect the PPT entangled state in Eq. ( 15.411) with Tr(W±p_) = — 2 ~ 2 V ^ — 
-0.29. 
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7.2 Approximated n-qubit Cluster SEWs 

For the odd case of cluster SEWs we add one of the statements ^-l^i ^° PH~6P as 

n" 

yap. ) 



k=0 



where m — 1, n". The eigenvalues of W'^ for all (ii,i2, v) e {0) 1}™ are 



(ap.) 



«o + + (-l) i2m a 2m + (-1)* 

j=0 



*2m+*2m+l , 



2m,2m + 



I ( 1 \«2m-l+«2m„ 

l T l ""J U 2m-1,2 



The coordinates of the vertices which arise from pure product vectors are listed in the Table 
5, where 



Product state 


(-Pi, ^3) •••) p2m-3) Plm-li P%m-, Pim+li P2m+3, ^W'+l, P2ro,2m+1> ^2m-l,2m) 


\{odd) 

h,*2. --.v'+i 


d>) 


((-l)* 1 , (-1)* 2 , -., (-l)^ 2 , (-l) im -\0, (-l)H (-l)V'+i,0,0) 


•l.^.-.V'+l 


$) 


, (_i)fa (-1)^,0, ±1, 0, (-1)^, (-l)W, 0, 0) 


p^ll{odd) 

n.*2.-.V'+l 


$) 


((_l)u (-1)^-2,0,0,0, (-1)^, .., (-l)V'+i, (-1)^,0) 


a ///(odd) 

»l,»2,-..,» n // +1 


|$) 


((_!)« o, 0, 0, .., (-!)*«"+! , 0, (— l)* m ) 



Table 5: The product vectors which seem to be the vertices for W 



,,(n) 



A /„(odd) = A ( orfd ) H (2m~2) M (2m-1) H (2m~l) M (2m) 

*1,»2. •••.»„//+! »l.»2.-.V/ + i 

By choosing any n" + 4 from 2 n " +1 + 3 x 2 n " points, the following half-spaces achieves 

\p2m + ( — l) n -P2m-l + ( — 1 Y' 2 P2m,2m+1 + {~^Y A P2m-l,2m\ < ^ 



(7.62) 



, m — 1 7^ = 0, n" 

for all (ii,i2j%) £ {0, 1}" (see appendix II). This shows that the problem does not lie in the 
realm of exact LP problems and we have to use approximate LP one. To do so, we shift 
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aforementioned hyper-planes parallel to themselves such that they reach to maximum value 
i] = -j=. On the other hand the maximum shifting is where the hyper-planes become tangent to 
convex region coming from pure product states and in this manner we will be able to encircle 
the feasible region by the half-spaces 



(7.63) 



s.t. 



\P21n + Plm-l + P2m,2m+1 + ^2m-l,2m| — ^ 
Pi < 1 

where again we have used the Clifford group and write just the generating half-spaces. Due to 
the above inequalities the problem is approximately reduced to 

minimize Tr(W'^ | t) (t I ) 

C (ap.) 

\P2m + ("l) n P 2m _l + (-l) l2 P 2m , 2m +l + (-l) i3 P2m-l,2m\ < 75 ,V (ii,i 2 ,i 3 ) £ {0,1}" 
\P2m + (-l) n P 2m +l + (-l)' 2 P 2m , 2m +l + {-l) h P2m-l,2m\ < ^ ,V (il,i 2 ,h) G {0,1}" 

|^ 2 fc+i|<l , k = 0,...,n" 

\P2m\ < 1 
|-P 2 m, 2 m+l| — 1 
|-P 2 m-l,2m| < 1 

(7.64) 

which can be solved by LP method. 

The intersections of the half-spaces in the above equation form a convex polygon whose vertices 
lie at any permutation of the coordinates -P 2m _i, P 2m , -P 2m , 2m +i, -P 2 m-i, 2m of the points listed in 
Table 12 and any permutation of the coordinates P 2m , -P 2 m+i, -P 2 m, 2 m+i, p2m-i,2m °f the points 
listed in Table 13 where P' 's are defined by shifting the P 's and for all (ii, i n , i 2 m, 2m +i) G 
{0,1}" +1 . 

Therefore, to be guaranteed the non- negativity of the expectation value of W'^ over all pure 



product states, the conditions 

% + E T io(-i) i2fc+1 « 2fc+1 + (-i) l2m « 2m + (-i) 



!2m,2m+l , 



2m,2m+l 



2— y/3 



(ap) 



a, , > 

2m — 1,2m — 



SUCh that (-l)^m-l + (_l)<2m _|_ (_ 1 )i2m-l, 2m + (_ 1 )i2m, 2m+ l = X 
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(pi pi pi pi pi pi pi pi pi pi \ 

l r H ^3' ^2m-3' r 2m-l' r 2m' r 2m+H ^2m+3' •"' ^2n"+l' r 2m,2m+li ^2m-l,2m/ 



((-l)* 1 , (-1)* 3 , (-l)* 2 " 1 - 3 , (-I)""" 1 , (-l) i2m , (-l) i2m + a , (-ljV'+l.f-l)^^!,^) 

9 ^2m + -^2m-l + ^2m,2m+l + ^2m-l,2m = ^ 



((-l)* 1 , (-l) i3 , (-l) i2m ~ 3 , (-1)* 2 " 1 - 1 , (-l)* 2m , (-l) i2m+1 , (-l) i2m + 3 , (-l) i2 «"+l, (-1)^.2-+!, ^=2) 

5 ^2m + ^2m-l + ^2m,2m+l + ^2m-l,2m = ~~ ^ 



Table 6: The coordinates of vertices for W'^ . 



(pi pi pi pi pi pi pi pi pi pi \ 

V-M) ^3' •■•> r 1m-%i r 1m-\i r 1mi r 2m+l ' r 2m+Z > *2n"+l> r 2m,2m+li r 2m-l,2m,l 



((-l)* 1 , (-1)^ , (-I)* 2 — 3 , (-1)^—1 , , (_l)<2m+l j (_l)i 2m +3 ; ..^ (_l)*2n''+l j (_ 1)*2 m ,2 m +l ? 2__V^ 

=) P' J- P' _|_ pi _|_ P' — 1 

^ r 2m ^ r 2m+\ ^ ^m^m+l ^ ^ 2771-1,2™ ~~ 1 



((-i)* 1 , (-i)S {-i) i2m - :i ,(-i) i2m - i ,(-iy 2m ,(-i) i2m + i ,(-i) i2m + :i ,..., (-iy^"+i, (-1)^,2^+1^ vj=2) 

3 ^2m + ^2m+l + ^2m,2m+l + ^2m-l,2m = ~~ * 



Table 7: The vertex points of approximated FR of W' 

°"(«P) 



; /(n) 

C7, 



a + Efclo(- 1 ) <2fc+lfl afc+ i + (-i)^^ + (-l)* 2m - 2m+1 a 2m , 2m+1 + ^a 2m _ lj2m > 

(7.65) 

such that (-l) i2 ™-i + (-I)m™ + + (_i)^ m , 2m +i = _j 

and any permutation of parameters a2 m -i, a2 m , a,2m,2m+i and a2 m -i, 2m together with the fol- 
lowing conditions 

% + E T io(-i) i2fc+1 «2 fc+1 + (-i) i2m a 2m + (-i) i2m < 2m+1 a 2m , 2m+1 + ^v-^ > o 

such that (-l) i2 — i + + (_i)<a».-i,2m + (_i)i2 m , 2m +i = x 

a + Efclo(- 1 ) i2fc+lfl afc+ i + (- 1 ) t2m « 2m + (-l) l2m ' 2m+1 a 2m , 2m+1 + ^^m-^m > . ccN 

(7.66) 

SUCh that (-l) i2 ™-l + (-1)«™ + (_l)Mm-l, 2m _|_ (_l)i2m,2m + l = _! 

and any permutation of parameters a 2m , a 2m+ i, a 2mi 2m+i and a 2m _i i2m . must be fulfilled. Sim- 
ilarly, one could repeat this approximate solve for the even case of cluster SEWs just like 
above. 
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Now we discuss the non-decomposability of cluster SEWs mentioned at the end of subsec- 
tion 15. 21 For more than three-qubit, the SEWs corresponding to half-spaces (I7.63P contain 
some nd-SEWs. As an instance, consider the SEW 

W = 4= he ~ S[ c) - SP - Sf^P - SPSP (7.67) 
V3 

The expectation value of W with the following density matrix 

P = i [i u + \{sp + spsp + spsp - spspsp)] 

is Tr(Wp) = —0.345 which means that W can detect p. On the other hand one can easily 
check that p is a bound entangled state and hence W is a nd-SEW. 

8 Conclusion 

We have considered the construction of EWs by using the stabilizer operators of some given 
multi-qubit states. It was shown that when the feasible region is a polygon or can be approxi- 
mated by a polygon, the problem is reduced to a LP one. For illustrating the method, several 
examples including GHZ, cluster, and exceptional states EWs have studied in details. The 
optimality and decomposability or non-decomposability of SEWs corresponding to boundary 
half-planes surrounding the feasible region have examined and it was shown that the optimality 
has a close connection with the common eigenvectors of SEWs. In each instance, it was shown 
that the feasible region is a polygon and the Hermitian operators corresponding to half-planes 
surrounding it are SEWs or positive. Also we have showed that, by using the Clifford group 
operations one can find vertex points and surrounding half-planes of feasible region only from 
a few ones. 
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Appendix I: 

Stabilizer theory 

Here we summarize the stabilizer formalism and its application to construct an interesting 
class of EWs so-called stabilizer entanglement witnesses (SEWs) |lUj . 

The / = 2 k (where k — 0, 1, ...,n) stabilizer states of n qubits can be thought of 

as representation of an abelian stabilizer group S n -k generated by n — k pairwise commuting 
Hermitian operators in the Pauli group Q n , which consists of tensor products of the identity 
I2 and the usual Pauli matrices a x ,a y and o z together with an overall phase ±1 or ±i [3, 
[25l [26] . The group iS n _fc has 2 n ~ k elements where among them we can choose S 1 , S n _ k as 
generators. This group leaves invariant any state in the stabilizer Hilbert space Hs spanned 
by {1^},-, |^)}, i-e., 

S\^) = \^) ,VSeS n - k ,V\i(>)eH s . (i-i) 

Similar to the Pauli matrices, for each element S of 5 n _fc the relation S 2 = I 2n holds and any 
two elements A and A. of this group satisfy 

Tr(AA) = / 2 „^. (1-2) 

The n-qubit Clifford group Clin) is the normalizer of Q n in U(2 n ), i.e., it is the group of 
unitary operators U satisfying UQ n U^ = Q n . It is a finite subgroup of U{2 n ) generated by the 
Hadamard transform H, the phase-shift gate M, (both applied to any single qubit) and the 
controlled- not gate CNOT which may be applied to any pair of qubits, 

/ 



1 1 
1 -1 





(1 o\ 


, M = 




J 





CN rs \j) r \k) s = \j) r \j + k mod 2} s . 
Generators of the Clifford group induce the following transformations on the Pauli matrices: 
H : a — > a , a — > —a , a 

x z ' y y 1 z 

M : a — > a., , a„ — >■ — o r , o. 



(7 , 



(1-3) 



a. 
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CW 12 : < 



I®<J X 
I®<J y - 
I <S> cr 



-(7,(8)(7. <r <g>7 



— cr (g) er 



(1-4) 



(1-5) 



and their actions on the eigenvectors of Pauli operators are 

H\x^) = |z±) 
M|x ± ) = 1^) 
Mi/|^) = |y±), 

In the following table, we give some examples of stabilizer groups together with the correspond- 
ing stabilized states (the states which are invariant under the action of the stabilizer group): 
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Stabilized state 


Generators of stabilizer group 


\Ao...o) 


S (GHZ) :=a W a m^ a (n) 

S (GHZ) . = a (k-i) a (k) fc = 2 ,3,..,n 




SP = oPaP 


|Cn> = ^r®2 = i(|0>. + |l>« 4 a+1) ) 


Sp = ai k ~%i%i k+1) fc = 2,3,...,n-l 
5(C) = ai™" 1 ^ 




o(Fi) _ ^(1)^(2) (3) (4) 
"J j ' — Ox O z O z <->x 


l*i Fi) ) = i£ se5 5|ooooo) 

Fi 


sp > = ^ 2 Vi 3 Vi 4 ^ 




5f ) = ^a^Vf 

5f ) = 




5 (Se) = ^) a f) a f)^) 










l*2 e) ) = ^E Weo l^) 

V 1 r / ^ ^ 






SW=aP*x 2) ...o-x 8) 


(*i»* 2 »* 3 ) G {Mi 3 


c(Ei) _ JIJJS)^) (5) (6) (7) 

o(Ei) _ J2) (3) (4) (5) (6) (7) 
"4 — Oz Oz Ox O x Uy Uy 

c(Ei) (1) (2) (4) (5) (6) (7) 

j\ ' = 0~ x 0~ x 0~ z 0~y 0~y 0~ z 






l,T f (Ni)\ _ / i /mnn\ _i_ n 1 1 ^ 

l*± ) •= (,T2^ 000 / ± l 111 ^ ) 


S<Ni) = (T W a W , S^)=afaf 


\<p) = ^ 2 (\v) + S 1 \v)) 


S (V) =S (GHZ) s ( v ,) =( g ) m iS (GHZ) 
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where 

£ ■= (10000000), |1010101), 10110011), 11101001), 10001111), 11100110), 11011010), 10111100)}, 
O : = {| 1111111), 11110000), |0100101), (1000011), |0010110), |0101010), 11001100), 10011001)}, 



X 2 = aPoPaPaP 



V - ^(1)^(4)^(5)^(6) 



') |0)i|0) 2 ...|0) 2m , ©& fc = 0, V (zi,i 2 ,...,i 2 J e {0,1} 

Here © is the sum module 2. 

Each of the stabilizer groups stated in the above table, corresponds to some graph states 
[36j [37] . These states are defined as follows: A graph is a set of n vertices and some edges 
connecting them. For every graph G, it is associated an adjacency matrix T whose entries are 
Tij = 1 if the vertices i and j are connected and T^- = otherwise. Based on that one can 
attach a stabilizer operator for every vertex i as follows 

The graph state \G) associated with the graph G is the unique n-qubit state satisfying 

S^\G) = \G), for i = l,...,n. 



2/ii 



The case k = and k > are called graph state and graph code respectively j26j [38]. One 
can denote the generators of any stabilizer group by a binary (n — k) x In stabilizer matrix 
\X\Z\ where X and Z are both (n — k)xn matrices. Matrices X and Z have a 1 whenever the 
generator has a and a z in the appropriate place respectively. For instance, in the five-qubit 
case, this form becomes 



/ 



[X\Z] 



10 10 
10 1 
10 10 
10 10 



110 0^ 

110 

1 1 

1 1 / 
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By adding k rows to the stabilizer matrix \X\Z\ such that the n resulting rows are linearly- 
independent, one can construct the matrix T (called generating matrix in coding theory) as 
follows 



X 


z 


X 


z 



\ 



For five-qubit in hand, we have 



r, = 











1 








1 


1 











1 








1 








1 


1 





1 





1 

















1 


1 





1 





1 





1 











1 

















1 


1 


1 


1 


1 



It is necessary to note that the added k rows are not unique and this freedom in choice leads 
to the several locally unitary equivalent graphs for a given graph code. 

By using the Gaussian elimination method on matrix T one can transform it to the standard 
form T' = [I\A], where I is a n x n identity matrix and A = T is adjacency matrix for the 
related graph G. The standard form of T 5 becomes 



/ 



1 
10 
10 
1 
1 

The graphs for the stabilizer groups stated in the above table are shown in Fig. 2. 

Appendix II: 

Proving the inequalities: 

For simplicity in the following proofs we introduce the abbreviations 



V 



1 1 
111 
10 1 
110 1 
11110 



\ 



/ 



a % := (ffW) A := (af) 7 , := (a 



Stabilizer EW 41 

+ # + 7- =1 (II-l) 

where (o^ ) with j = x,y,z are the expectation values of Pauli operators on any arbitrary 
pure qubit state. 
The proof of (SHI]): 

We give the proof only for the j = 2 since the proof for other cases is similar. 

n> n' 

\ p i ±p 2 + ID-i^j < i^i + 1^1 + E ^ 

fc=l fc=l 



l«l« 2 l I l«3-«„l + E l a 3«4---« 2fc -2A,2 fe -lA,2 fc « 2fc+ l---«nlJ + I P1P2 ^3 -«„ I + ITi7 2 | < 

l«i« 2 l + \PiP a <x a ...a n \ + IT1T2 1 < l«i« 2 l + + I7i7 a | < 1 

The last inequality follows from the Cauchy-Schwartz inequality and (IH-ip . 

To show that the big bracket is smaller than one, we write 

w 

^•••^nl ^ \ l a 3 Q; 4--- a 2fc-2'^l,2fc-l/^l, 2 fc a 2 
k=2 

= \a 3 ...a n \ + \(3 3 p 4 a 5 ...a n \ + \a 3 a 4 p 5 f3 6 a r ...a n \ + \a 3 a 4 a 5 a 6 (J 7 (3 s a 9 ...a n \ 

+ ... + |a 3 a 4 a B ...a an ,_ a ^ aB ,_ 1 /3 an/ | < \P 3 p 4 a 5 ...a n \ + \a 3 a 4 \ 

x [|a B ...a n | + \P B P B a 7 ...a n \ + \a 5 a 6 (3 7 (3 s a 9 ...a n \ + ... + \<x B a e ...a 2nl _ a P ;ia ,_ 1 P 2a ,\] 

< \P 3 P 4 a B ...a n \ + \a 3 a 4 \[\(3 5 (3 6 a 7 ...a n \ 

+|a 5 a 6 |[|a 7 a 8 |[...[|^ 2n ,_ 3 /3 2n ,_ 2 a 2n ,_ 1 a 2n ,| + W 2nl _ 3 a 2n ,J[\a 2nl ^a 2nl | + \P 3n , _&„,]]...]]] 

By the Cauchy-Schwartz inequality, the last bracket [|ol , ,a: ; | + \P ,_,P„ ,] is less than or 
equal to 1. Replacing it with its maximum value 1, the same argument holds for the term 

[l/ 3 2 n'-3^2n'-2 a 2n'-i a 2n' I + I "w-s Vj 1 1 ■ Proceedin g in tnis wa Y, w e deduce that the big bracket 
is less than or equal to 1. 
The proof of flX24ft : 

l7 2 m-2 < -^2m-l7 2 m T 2 m-3 ^2m-2 72m- 1 7 2 m-2 Am- 1 Am 72m + l I — 
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lT2m-2 < ^2m-lT2m I I I'im-Z ^2m~2 12m- 1 I ~l~ I T2m-2 Am- 1 Am Tsm+l I 

taking 7 2m „ 3 = 7am+1 = 1 

— I T2m~2 ^2771-1 T2m I I ^2771-2 O^m- 1 I I T2m-2 Am- 1 Am I 

Cauchy-Schwartz inequality yields 

< ( a 2 + 2 +/5 2 \*/ 2 2 + 2 + 2 o2 \ § 

— V 2m-l '2m— 1 "2m-U V 1 2m-2 ' 2m 2m-2 '2m-2"2m/ 

< f 7 2 (J +fl 2 )+o J < 1 

— V '2m-2V '2m ~ ^ 2m J ~ " 2m _ 2/ / — - 1 - 

In the above inequalities we have used and the fact that 



aa? + bft + cry? < 1 (II-2) 

when the positive coefficients a, b, c are less than or equal to one. The proof for even case is 

the same as the odd case. 

The proof of ffFSTb : 

We prove only the case 1 = 1. 

n' 

k=l 

= | cos(0j_i) cos(9j) + sin(6 l 1 ) sin(# 2 )... sin(6* n ){cos(0i) cos(0 2 )... cos(<^ n ) 

n' 

+ ^ cos (0i) cos(0 2 )--- cos(0 2 z c _ 2 ) sin(0 2fc _i) sin(0 2 fc) cos(0 2fc+ i)... cos(0 n ) 

k=l 

+ cos(0i) sin(0 2 ) sin(0 3 ) cos(0 4 )--- cos(</> n )}| 
< | cos(^_i) cos(#j)| + | sin(^_i) sin(^)| x | cos(0i) cos(0 2 )... cos(<p n ) 

n' 

+ ^ cos (0i) cos(0 2 )... cos(0 2fc _ 2 ) sin(0 2fc _i) sin(</> 2A .) cos(0 2fc+ i)... cos(0 n ) 

k=l 

+ cos(0i) sin(0 2 ) sin(0 3 ) cos(0 4 )... cos(</> ra )| 
On the other hand, we note that 

n 1 

cosOi) cos(0 2 )... cos(0 n )+^ cos(0i) cos(0 2 )- cos(0 2fe _ 2 ) sin(0 2jt _i) sin(0 2fe ) cos(0 2fc+ i)... cos(0 n ) 



fc=i 
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1 ~h 

+ cos(0i) sin(0 2 ) sin(0 3 ) cos(0 4 )...cos(0„) < — - — 

Hence we get 

IP, + P x + £ P ^ + P»\ < I cos(^-i) 008(^)1 + MOj-i) sm(%)| < 

k=l 

The proof of < TT62l} : 

The proof is for m > 2 and the half-spaces with positive coefficients. The proofs for the other 
cases are similar. 

\P2m + Plm-l + P 2 m,2m+1 + P2m-l,2m| = 

lT2m-l ( -^2mT2m + l T2771 -2 ^2m- 1 T2771 "I" T2m-1 /^2m/^2m+l T2m+2 T2m-2 l^2m- 1 Am T2m + 1 I — 

I cos(# 2m -i) sin(6> 2m ) cos(0 2m ) cos(6» 2m +i) + cos(# 2m _ 2 ) sin(6> 2m _i) cos(0 2m -i) cos(0 2m ) 

+ cos(# 2m _ 2 ) sin(6> 2m _i) sin(</> 2m _i) sin(6> 2m ) sin(0 2w ) cos(0 2m+ i) 

+ cos(6> 2m _i) sin(6> 2m ) sin(0 2m ) sin(6> 2m+1 ) sin(0 2 m+i) cos(6> 2m+2 )| 

We note that the maximum value of the right-hand side of the above statement is ^= . Hence 
we get 

|P2m + P2m-1 + P 2 m,2m+1 + P2m-l,2m| — ~7= 

Appendix III: 

Odd case of GHZ SEWs 

Let us consider the Hermitian operator 

n n" 

W'H = a L n + y a k S< GHZ > + V a 2L+1 S^S^ (III-l) 

GHZ 2" / j k k / j 1,2/c+l 1 2fe+l v ' 

fc=l fe=l 

coming from (13.71) by replacing all even terms with odd ones S[ GRZ > Sf^*^ 2 ' (the name odd refer 
to the index 2k + 1). It is easily seen that the eigenvalues of W are 

n n" 

«o+E(- 1 ) iia i+E(- 1 ) il+i2fe+la ^ +1 > v(M 2 ,...,ye{o,ir (ni-2) 
i=i fc=i 
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(Pi Pi P 1 P 


P 1) -fl,3j -Pi, 5) • 


• • , P\,2n"— 


1 ! 1 l,2n"+l / 






(0,0,... 


,0,0, ±1,0,0, 


...,0,0) 








(0,0,... 


,0,0,0, ±1,0, 


...,0,0) 




A' |tf±) 




(0,0,... 


,0,0,0,0,0,... 


,0,±1) 




Hi 2 ,...,iJ* + ) 


((" 


l)<a,(-l)V* 8j .„,(- 


l)*„-2 


-1 )*„-!+*« 


,0,0,0,. ..,0,0) 



Table 8: The product vectors and coordinates of vertices for W'^ z . 
The product vectors and the vertex points of the feasible region are listed in the table 8, where 

A' fc = (M^) j M^ k+l \ k = 1,2, n" 

and \§! ) is defined as in (13.101) . The convex hull of the above vertices, the feasible region, 
comes from the intersection of the half-spaces 

n" 

\P 1 ±P. + J2(-^) ik P i ,, k J<^ , j = 2,...,n , V(h,i2,...,in»)e{0,n n " 



(III-3) 



k=l 



and is a [n — l)2 n " +2 -simplex. 

We give the proof only for the j = 2 since the proof for other cases is similar. 



IP, ±p 2 + J2(-i) lk P iak+1 \ < \pa + \p t \ ±Ei p ^i 



< 



k=l 



k=l 



K« 2 I K-aJ+X)! 

tt 3 tt 4 .. ■Oi 2k _ 1 



k=2 



|«i«al + IAA« 3 --«J + |7i7al ^ K« 2 I + + I7i7 2 l < 1 



The last inequality follows from the Cauchy-Schwartz inequality and (1II-lj) . 
Now the problem is reduced to the following LP problem 



minimize = a + ££ =1 a k P k + ELi a li2fc+1 Pi, 2 fe+i 

subject to \P 1 ±P j +J2£i(-l) ik Pi,w\<l, j = 2,...,n , V (n, i 2 , vO 6 {0, 1}"" 

(III-4) 
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By putting the coordinates of vertices (see table 2) in the objective function J- W ,( n ) and 



GHZ 



requiring its non-negativity on all vertices, we get the conditions 

% > , a > |aj , a >E"=2kl 



(III-5) 



a >k, 2fe+ il k = l,...,n" 



on parameters <2j. 

Odd case of cluster SEWs 

Let us consider the Hermitian operators 



< n) = «o^+E^ +1 C 1 +«-^ ) + ^ +1 ^ ) ^ 1 > m = l,..,n» (IH-6) 



fc=0 



Note that instead of the last term we can put the term a, , „ Sjp^ with m — 1, ...,n'. 

Jr 2m— 1,2m 2m— 1 2m ' ' 

The eigenvalues of are 



«« + E(- 1 ) i2j+1 ^ +1 + (- 1 ) i2m °- + 



«2m+*2m+l . 



l 2m,2m + 1 , V (* X , 2 2 , . . . , 2 n ) G { , 1 }" 

(III-7) 

The product vectors and the vertex points of the feasible region are listed in the table 9 



Product state 


(-Pi) -f3) •••) P2m-3) P2m-1) -FWi) ^m+l) -P2m+3) •••) ^W'+lj -P2m,2m+l) 


^(odd) |^\ 

»l,«2.-.«„/'+l 


((_!)*, (_1)« (_l)^- 2 , (-l)^-S 0, (— l) im , (-l)*™+\ (-l)V'+i, 0) 


A /(odd) 

il,»2, -,<„//+! 1 ' 


((_l)u 5 (-1)^,0, ±1, 0, (-l)Wi, (-l)W, 0) 


^//(odd) iq\ 

»1,»2. -.*„"+! ' 


((_!)* (_i)*2 > (-1)^-2, 0, 0, 0, (-l)Wi 5 (_l)V' + i, (-1)^) 



Table 9: The product vectors and coordinates of vertices for W 



/(n) 
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where 

|$) = |^) 1 |x+) 2 |^) 3 |x+) 4 |,+) 5 ...|x+) n _ 1 |z+) n 

A !;t 3 ,„ = ®u hU) > v(i 1I < 2l ...,w)G{o,ir+ i 

j^i(odd) _ j^(odd) o-(2m-l) rr(2m) o-(2m+l) 

ii,»a.-.* n /'+i 'i'^' ■•••V'+l 

^//(odd) = ^(odd) H (2m-1) M (2m) H (2m) M (2m+1) 

<l.*2.-.*„/'+l *1»*2 V'+l 

For a given m, the feasible region (the convex hull of the above vertices), comes from the 
intersection of the half-spaces 

\P2rr1 i ^2m-l + -P2m,2m+l| < 1 
|^2m =t ^2m-l ~~ ^2m,2m+l| — 1 
l-f^m =t ^2m+l + -P2m,2m+l| < 1 

(III-8) 

\P21n i -P2m+1 — P2m,2m+l\ — 1 

|-P2fc+i| < 1 , m,m-l^fc = l,...,n" 
and is a (2n" + 12)-simplex. For the first two inequalities we have 

ll^m— 1 ^2ra72m+l T2m-2 ^2m~ 1 T2m ~l~ T2m— 1 Am^2m+l72m+2 I — 
l72m-l ^2m72in+l I \l 2m-2^2m-l1 2m\ 1 7 2 m- 1 Am (^2m+l 1 2m+2 I 

taking 7am _ 2 = 7 2m+2 = 1 

— I Tiro- 1 ^2t72m. + l I 4~ l < - t 2m-lT / 2m I 4~ \l 2m-\ (^2m(^2m + l\ — 

Cauchy-Schwartz inequality yields 

< ( a 2 +7 2 + «2 / 2 2 + 2 + 2 ^2 \£ 

— V 2m '2m "2m/ V <2m-l '2m + l 2m-l '2in-l"2ni+l/ 

<(7L- 1 (7L + 1 +/3L +1 )+«2m- 1 )"<l 

In the above inequalities we have used ( 111-11) and ( 111-21) . 
For the second two inequalities we have 

lT 2 m-l < -'-2mO'2m + l 72m ^2m+l 72m+2 72m- 1 Am /^2m+l 72m + 2 I — 
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lT2m-1^2mT2m + l I > I T2m ^2m+l T2m + 2 I ~T~ I T2m- 1 H2m P^m + l T2m + 2 I 

taking 7am _ 1 = 7am+2 = 1 

— I^SmTam+l I \l ImP^lm+lA |/^2m/^2m + ll — ^ 

The last inequality follows from the Cauchy-Schwartz inequality. 
The objective function is 

n" 

Fyvg) =a + Yl a »+i P 2fe+i + a ^ P ^ + a 2™,2 m+ i P 2 m ,2 m+ i , m = l,...,n" (III-9) 

where 

P 2k+i =Tr(S£l Ps ) , P 2m , 2m+1 =Tr(SWsW lPs ), 

If we put the coordinates of vertices (see table 4) in the objective function (1III-9I) and 
require the non-negativity of the objective function on all vertices, we get the conditions 

> S^ n " I I 

a — 2-uj=0 \ a 2j + l I 

a o — Sj=l l a 2j+il + Y2j=m+1 \<hi+i\ + kJ (iii-io) 

a — 2—ij=\ l fl 2j + ll ' Z^j = m+1 l a 2j + ll ' l fl 2m,2m + ll 

for parameters a.j. 
Exceptional SEWs 

Here we mention briefly the SEWs that can be constructed by the stabilizer operations of 
the five-qubit, seven-qubit, eight-qubit, and nine-qubit states that can be solved by exact LP 
method. 

Five-qubit SEWs 

Consider the following Hermitian operator 

W FI = a I 25 + ai S^ + a 2 SP + a 3 SP + a 3>4 S^S^ 
Eigenvalues of W F4 are 

3 

a + ^(-l)H -±a 3 ,4 , V (ix, is) E {0, l} 3 

3=1 
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Product state 


(Pi, P2, P3, ^3,4) 


H^H^ (SW) 15 |^f 1 )) 
P( 3 )p( 4 ) (SW) 25 |^f 1 )) 
HUH® (SW)JVW) 


(±1,0,0,0) 
(0, ±1,0,0) 
(0,0, ±1,0) 
(0,0,0,±1) 



Table 10: The product vectors and coordinates of vertices for W Fj . 

The vertex points of the feasible region are listed in Table 10 
where 

\9™) = \^) 1 \*+) a \*+) a \x+) < \ > B 
(SW) y = (CNUCNUCN)^ 
The operator (SW)ij when acts on any two arbitrary pure states swaps them, i.e., (SW)ij\^)i\ip)j 
\ip)i\ijj)j. Inequalities obtained from putting the vertex points are 

o > |oJ % = 1,2,3 , o > |a 34 | 

Boundary half-spaces of feasible region are 



(III-ll) 



\P ± P + P + P I < 1 IF ± P + P - P I < 1 

|- 1 l- 1 -- 1 2^- I 3^ 3,4l — x 1 \ 1 1 - 1 - 1 2 ~ 1 3 3,4 I — 

|P±F-P+P I < 1 IP ± P - P - P I < 1 

I - 1 1 — J 2 x 3 ~ 3,4 I — ) |-'l- 1 -- I 2 J 3 3,4l — 

We prove only the following inequality since the proof of the other inequalities is similar to 
this one. 

\ p i + P2 + P3+ P 3 a\ = K727 3 a 4 + «27 3 74«5 + aia 3 747 5 + Pia 2 a 3 f3 A < 
l a i7 2 7 3 Q! 4 | + K73 7 4 a 5 l + Ka 3 74 7 5 l + \Pia 2 a 3 (3 4 \ < 



I«i72« 4 | + |a2 74a 5 l) + l« 3 KI«iT 4 T 5 1 + \Pia 2 p 4 \) < 
\% I K + ll) 1 (ay, + 7 X) * + \a s I (a? + #) 1 (ihl + a 2 ft) * < 

l7 3 l + 7X) 1 + l«sl (7 4 2 7' + ^AY < 



Stabilizer EW 49 



(«s + 7a ) 1 (*X + 7X + 7 4 2 7 5 2 + « 2 2 /? 4 2 ) * < 

{« + 7 4 2 (« 2 + 7 2 ) + 1 < (aX + 7 2 + a 2 /? 2 ) * < (a 2 + 7 2 + /3 2 ) * < 1 

The above inequalities follow from the Cauchy-Schwartz inequality and the equations (1II-lj) 
and (ED - 

From 2 4 eigenvalues of W Fi , six of them can take negative values. If we take all a 1 ,a 2 ,a 3 ,a 34 
positive and without loss of generality assume that a 1 > a 2 > a 3 > a 34 , then these eigenvalues 
are 

a — a 1 — a 2 ± a 3 + a 3 4 , a — a 4 — a 2 ± a 3 — a 3 4 , a ± a 1 =p a 2 — a 3 — a 3 4 

Seven-qubit SEWs 

Consider the following Hermitian operator 



« i 1,4 1 4 

In addition to the above operator, we can consider other Hermitian operators which differ from 
the above operator only in the last term, that is the last term of them is one of the following 
operators 



1 4 ' 2 5 ' 3 6'1 5 ' 2 



Eigenvalues of W Se are 

6 

% + £(-i)S + (- 1 ) il+ ' 4 %,4 v (*i> ->^) e {0, l} 6 

5=1 

The vertex points of feasible region are listed in table 11 
where 

\<S>(S*)) = \z+) 1 \z+) 2 ...\z+) 7 
A ffi^ = (^ 1 ^( ff * ) )* a ( ff * 4) )*" 

A tL = (^ 1) ) <1 (^ ) ) <a (^ ) ) <> ®;=i^ ) . v ( il ,i 2> < 8 ) e {o,i} 3 



Af°) = (ai 1 ^ ®J =1 MW-VHW-V, V 2 e {0, 1} 
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Product state 




A(Se) |<J>(Se)\ 




A'( Se ) |$( Sc )\ 


(0, 0, o, (-i) fa , (-1)* 8 , o) 


A(So)|$(Se)^ 


(0,0, 0,0, 0,0, (-1)*) 



Table 11: The product vectors and coordinates of vertices for W Se . 
Boundary half-spaces of feasible region are 

l^iP. + Pj <1 , \R±R-P 1A \<1 z = l,2,3 , j = 4,5,6 (HI-12) 

Although all of the inequalities (1III-12I) can be derived by Cauchy-Schwartz inequality but 
require a tricky way. The proof of two cases % = 2,j = 6 and i = 3, j = 5 are similar and 
therefore we prove only the former case. 

|7 2 73 7 6 7 7 I + \a 4 a 5 a 6 a 7 \ + |&/3 3 /? 5 /3 7 | 

taking 7 2 = a 4 = j3 x = 1 

< l7 3 7 6 7 7 l + |a 5 a 6 «7l + WsPrl < 

bWMl + ft + ft) + ftft(ft + ft + ft) + filial + ft + ft)} * < 

[(ft + ft + ft) (ft + ft + ft) (ft + ft + ft)Y=l 

The above inequalities follow from the Cauchy-Schwartz inequality and the equations (1II-lj) 
and (ED - 

Inequalities obtained from putting of vertex points are 

a > |a M | , a > Y? j=1 |aj , a > E;=4 Kl 
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Eight-qubit SEWs 

Consider the following Hermitian operator 



Eigenvalues of W Ei are 



^+E(- 1 )^, + (- 1 ) n+l2+J3a M,3 
J'=l 



Ul+i2+*4, 



"1,2,4 



V (ii,i 2 ,...,i 5 ) e {0, l} 1 



The vertex points of feasible region are listed in table 12 



Product state 


(P P P P P P P ) 

V- 1 11 ± 2 1 x 3 > 4.1 5 1 1,2,3' 1,2,4/ 


|$(Ei)> 


(±1,0,0,0,0,0,0) 


HWHW...HW\$& r )) 


(0, ±1,0, 0,0, 0,0) 


hWmWh®mw\$<&)) 


(0,0, ±1,0, 0,0,0) 


hWhWmWmW\$<&)) 


(0,0,0, ±1,0, 0,0) 




(0,0, 0,0, ±1,0,0) 




(0,0, 0,0,0, ±1,0) 




(0,0,0,0,0,0,±1) 



Table 12: The product vectors and coordinates of vertices for W Ei . 



where 



|$(Ei)\ 



\x ± ) 1 \x + ) 2 ...\x^ 



Choosing any seven points among the above vertices give the boundary half-spaces surrounding 
the feasible region as follows 



\p, + (-i)h p 2 + (-1)^3 + (-iy*p 4 + (-i) i4 ^ 5 + (-l) 45 ^,: 

, V(2 1 ,Z 2 ,...,2 6 )G{0,1} 6 



-iy 6 p i < i 



(111-13) 
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We prove only the following inequality since the proof of the other inequalities of ( 1III-13I) is 
similar to this one. 

\p + pj r p + p + p + p + p | = 

|- I l^- I 2^- t 3^- , 4^- t 5^- , 1,2,3 ~ ± 1,2,4 I 

|a 1 a 2 ...a 8 +7 1 7 2 ...7 8 +7 1 a 2 /3 3 7 5 7 6 /? 7 +7 2 7 3 a 4 a 5 /? 6 /? 7 +^^ 

l«2l(l«il(« 4 2 + J 2 A )Ha 2 3 a 2 s ...a 2 8 + P 5 2 /3 6 2 ^ + \%(3 3 %%P7\ + |/3ia 3 7 4 7 5 /38|) + 
\% I (« 4 2 + P. 2 + ll) 1 {lllH-ll + 7 2 « 2 /V/V + afa^A. 8 ) 1 < 
l« 2 1 (« 2 + /?i 2 + I 2 ) 1 (« 2 « 2 -« 2 + /W^ + /?3 2 7 5 2 7 2 /5 7 2 + « 2 74 2 7 2 /? 8 2 ) h + 
l7 2 l (T^Tf-T^ + 7 3 2 ^/?6 2 /? 7 2 + ^^^/? 8 2 )^ < (al + 1 lY x 
(a 2 a 2 ...a 8 2 + /? 5 2 /3 6 2 7 7 2 + /?3 2 7 5 2 7 2 /V + + 7 1 2 7 2 7 5 2 ...7 8 2 + lK^^ + 1 < 

[oj (« 6 2 (« 2 « 2 «s + <W + 7 2 /?6 2 /? 7 2 ) + d, 2 ^ 2 ! 2 + 7 5 2 (/?3 2 7 2 /3r 2 + « 2 7 4 2 /V + 7 2 7 2 7 2 7 2 7 2 )] 1 < 

K + /5 5 2 + 7 2 ) l <i 

where, we have used the Cauchy-Schwartz inequality and the equations (1II-lj) and (1II-2j) re- 
peatedly. 

Inequalities obtained from putting the vertex points are 

% > |aj i=l,...,5 , a >|a 1>2i3 | , a > |a 1A4 | 
Nine-qubit SEWs 

Consider the following Hermitian operator 



Eigenvalues of W Ni are 



l o + + (-l) n+l3 ai,a V (« ls i 2 , i 8 ) G {0, 1}' 
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Product state 


/ p pppppppp\ 

K 1 1 ) J 25 3' 4' 5' 6 ' 7' 8 ' 1,3/ 


A (Ni) |«f,(Ni)\ 




A (Ni) |flj(Ni)\ 

»li*2i*3 


U~ -U 1 U, U, U, U, U, ( — Ij ,(-lJ J ,UJ 


A/(Ni) |^,(Ni)\ 


l u > l — L ) A~ L ) >\~ L ) 5 u, u, u, U, UJ 


A //(Ni) |$(Ni)\ 

n.i2.*3 


(0,0,0,0,0,0,(-l)M-l)M-l)^) 


A (Ni) i $ (Ni)\ 

i l. i 2> i 3' i 4. i 5' i 6 


(0, 0, (-rr, (-1)* 2 , 0) 



Table 13: The product vectors and coordinates of vertices for W Ni . 

The vertex points of feasible region are listed in table 13 
where 

|$(N0) = |x+) 1 |x+) 2 ...|x + ) 9 

Aj2 = (rf ) ) <1 (^) < ' 
A(Ni) , s = (aPY^aPY^a^H^H^H^ 

A'^;) , 3 = (ai 1) ) h (^ 3) ) i2 (^ 4) ) l3 ^ (1) ^ (2) ^ (3) 

A"^). 3 = ( ( 7i 1) ) il ( ( 7i 7) ) i2 ( ( 7i ; 9) ) i3 (M( 1 ))tM( 2 ) 

A 22^ = (^) <2 (^) <3 (^) i8 (^) i6 ®U ^ 

which in all of the above operators we assume that G {0,1} : ', with j = 2,3,6 . 

By choosing any eight points among the above vertices give the half-spaces surrounding the 
feasible region as follows 

<1 , l^-PiPj^l , 2 = 3,4,5,6 
\P 2 ±R\<1 , j = 7,8 
The proof of the above inequalities are straight forward. Inequalities obtained from putting 
the vertex points are 

a > + l a 2 l ) a o — l a J + l a 7l + l a sl 
a > |a 2 | + |a 3 | + |a 4 | , a > Xj= 3 l a J > a > |a 7 | + |a 8 | + |a 13 | 
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Figure Captions 

Figure-1: 8-simplex displaying the feasible region of the two-qubit GHZ SEW. 

Figure-2: Graphs corresponding to different graph states where the first two ones are 
graph states and the others are graph codes, (a) The star graph describing a GHZ state, (b) 
The linear graph describing a cluster state. The graph codes for (c) five-qubit , (d) seven-qubit, 
(e) eight-qubit and (f) nine-qubit stabilizer groups. 



